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Abstract 



Ph We prove a central limit theorem for an additive functional of the d-dimensional frac- 

tional Brownian motion with Hurst index H G (5+3) 3)) using the method of moments, 
extending the result by Papanicolaou, Stroock and Varadhan in the case of the stan- 
dard Brownian motion. 

Keywords: fractional Brownian motion, central limit theorem, local time, method of 
moments. 

k*" ■ Subject Classification: Primary 60F05; Secondary 60G22. 

vo 

CN 

^°. : 1 Introduction 

o 
en 



Let \B(t) = (5 1 (t), . . . , B d {t)),t > 0} be a d- dimensional fractional Brownian motion (fBm) 
with Hurst index H e (0, 1). The local time of B, denned as L t (x) = J 5(B(s) — x)ds, 
for t > and x e M. d , where 5 is the Dirac delta function, exists and is jointly continuous 
in t and x if Hd < 1 (see [2]). For any integrable function / : M. d — > R., using the scaling 
property of the fBm and the continuity of the local time, one can easily show the following 
convergence in law in the space C([0, 00)), as n tends to infinity 

(n Hd ~ l J™ f(B(s))ds,t>0) 4 (L t (0)Jj(x)dx,t>0y (1.1) 

If we assume that J Rd f(x) dx = 0, a central limit theorem holds with a random variance. 
In order to formulate this theorem, we need to introduce some notation. Fix a number 
f3 > and denote 

H{j = If e L\R d ) : f |/(x)||x| /3 rfx<cx) and / f(x)dx = o\. 
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For any / G H(f, and assuming (3 G (0, 2), the quantity (see Lemma 4.1 in [3]) 

11/11?:=-/ f(x)f(y)\x-yfdxdy = c M [ \f(x)\ 2 \x\-^ d dx (1.2) 

JR 2d JR d 

is finite and nonnegative, where / denotes the Fourier transform of /. 
Then, the following central limit theorem holds. 

— — d 

Theorem 1.1 Suppose -^ < H < ~ and f G H" . Then 

£ f(B(s))ds,t>0) ^ ( V ^r d \\f\\ JT _ d W(L t (0)),t>0) 



n 2 



in the space C([0,oo)) ; as n tends to infinity, where W is a real-valued standard Brownian 
motion independent of B and 

This theorem has been proved by Hu, Nualart and Xu in the reference [3], in the case 
where the Hurst parameter H satisfies -t-j- < H < -j, and it has been conjectured in that 
paper that the result can be extended to the case -^ < H < -^. The purpose of the 
present paper is to prove this conjecture. With this aim we will develop a new approach to 
prove Theorem 1.1 based on Fourier analysis. 

Note that the lower bound j-^ is optimal because for H < -^ the constant Cn,d is 
infinite. When d — 1 and H — |, the above theorem was obtained by Papanicolaou, 
Stroock and Varadhan in [41 with Ci, =2. 

2' 

As in the reference [3], the proof of Theorem 1.1 is based on the method of moments. 
In order to handle the integrals on [0,t] 2m , with respect to the measure dsi ■ ■ ■ds2 m , we 
make the change of variables «2fc-i — ^( s 2fc — S2k-i) and u^k = S2k, I < k < m. Then, the 
increments of B in small intervals will be responsible for the independent noise appearing in 
the limit. The main novelty of our approach, in comparison with [3], is a new methodology 
based on Fourier analysis and an iterative procedure in order to get the right estimates 
to derive the tightness of the laws and to show the convergence to zero in the truncation 
argument. 

After some preliminaries in Section 2, in Section 3 we prove some technical estimates 
based on Fourier analysis which play a fundamental role in our approach. Finally, Section 4 
is devoted to the proof of Theorem 1.1. Throughout this paper, if not mentioned otherwise, 
the letter c, with or without a subscript, denotes a generic positive finite constant whose 
exact value is independent of n and may change from line to line. 



2 Preliminaries 

Let {B(t) = (B l (t), . . . , B d (t)),t > 0} be a <i-dimensional fractional Brownian motion with 
Hurst index H G (0, 1), defined on some probability space (Q, J 7 , P). That is, the compo- 



nents of B are independent centered Gaussian processes with covariance 

E {B\t)B\s)) = ^(t 2H + s 2H -\t- s\ 2H ). 

The next lemma (see Lemma 2.1 in [3]) gives a formula for the moments of the increments 
of the process {W(L t (0)),t > 0} on disjoint intervals, where W is a real-valued standard 
Brownian motion independent of B. 

Lemma 2.1 Fix a finite number of disjoint intervals (a$, h] in [0, oo), where i = 1,...,N 
and bi < Oj+i. Consider a multi-index m = (mi, . . . , rajv), where m>i > 1 and 1 < % < N. 
Then 

N 

E(l[[W(L b M)-W(L ai (0))P) (2.1) 

t=i 

( II — 1& ) L m . det(A(w))~* dw if all mi are even 

Vj=l 2T- (27r)T-(T7ii/2)!/ ./nKA]^ 1 

»=i 

otherwise, 

where A(w) is the covariance matrix of the Gaussian random vector 

(B(w l k ) :\<i<Nandl<k< — -J. 

As a consequence, the law of the random vector (^(1/^(0)) — W(L ai (0)) : 1 < i < N) is 
determined by the moments computed in the above lemma. 

We shall use the following local nondeterminism property of the fractional Brownian 
motion (see [1]): For any n > 2 there exists a positive constant ku depending on n, such 
that for any = so < Si < ■ ■ ■ < s n < oo and U\, . . . , u n e M. d , 

n n 

Var ( Yl u i ■ ( B ( s i) ~ B ^i))) >knJ2 H^ ~ S ^ H " ( 2 ' 2 ) 

3 Technical estimates 

We are interested in the sequence of stochastic processes defined by 

fat 

F n (t)=n E ^ 1 / f(B(s))ds. 
Jo 

For < a < b < oo and m e N, let /" = E [(F n (b) - F n (a) m }. It is easy to see that 

„ m. 

r m = m \n miM ^ I E (lJ/(£(*))) & 

„ „ m -. m 

/ / ( n /(^)) ^ ( - 2 var ( e ^ • 5 ^))) ds ^' 
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where c m4 = j^a and D m = {(s t , ..., s m ) : na < s ± < ■ ■ - < s m < nb}. Making the 

m 

change of variables Xi = Y2 Uj (with the convention that x m+ i = 0) we can write 



c m ,d n 2 



/ / ( ] I f( Xi ~ Xi +^ 

.. m 

exp ( - -Var (^g r (-B(sj) - B(si-i)))j dsdx. 



The main idea in order to estimate these terms is to replace each product f($2i-i~ %2i)f{%2i~ 
%2i+i) by f{-X2i)f(x 2i ) = \f(x 2i )\ 2 . Then, the differences f(x 2i -i - x 2i ) - f{-x 2i ) and 
f(x 2i —x 2i+ i) — f(x 2 i) are bounded by constant multiples of |x 2 j_i| Q and |x2i+i| a , respectively, 
for any < a < 1, because /(0) = due to the fact that / has zero integral. We are going 
to make these substitutions recursively. To do this, we introduce the following notation. 



Let I^o = Im- For k = 1, . . . , m, we define 

JK™ d J Dm .•_!. , 1 



■2-i ^j+lj 



i=fc+l 



t=l 

where 



exp ( - -Var (^j r (-B(sj) - £(sj_i)))j rfs dx, 



/ 
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II 1/0*^)1 /(-Zk+i), if fc is odd; 



2 
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I! 1/(^)1 , if k is even. 

The following proposition controls the difference between I^ lk _ 1 and i^ fc . We fix a positive 
constant 7 such that 

.'±=f^ if I -Hd<H; 

\2H-1-Hd iiH<1 - Hd<2 H. 

Proposition 3.1 For fc = 1,2, . . . ,m, there exists a positive constant c, which depends on 
7, such that 

\jn _ jn I < rrl -J(U _ n \ ml =T^-l 

\ 1 m ) k-l 1 m,k\— Cn \° a ) 2 

Proof. The proof will be done in several steps. 

Step 1. Suppose first that k — 1. Applying the local nondeterminism property (2.2) and 
making the change of variable u\ = sx, Ui = Si — Sj_i, for 2 < % < m, we can show that 
I-PJjo — Imx\ i s l ess than a constant multiple of 



„ „ m m 

S **~ k / \f(xx -x 2 ) -f(-x 2 )\( 1 T|/(xi-Xj+i)|) exp ( - -£- y^\xi\ 2 u\ H ) dudx, 



4 



where 

O m = {(ui,. ■ ■ ,u m ) : < Ui < n(b- a), i = 1, • •• ,m}. 

Taking into account that that |/(x)| < c a |x| Q for a G [0, 1], we obtain 

„ „ 771— 1 

I rn _ jn l < •n mS ^T L I I 1-r, \ ai 1 T (\r-\° li 4- l-r- J^MIt l Qm 

\ 1 m,0 J m,ll^ c l n 2 / / \ x l\ \_[\ X i\ + F't+ll J Fm| 

771 

>, F»l M 7 ) " M " x 

7 = 1 

„ „ m— 1 

= Cl n m ^W / N ai ]I(l^r^mr0k7nr (3.2) 

s jR md JO m i=2 

771 

x exp f y. \ x i\ 2 ' u 'i H ) °^ dx, 

7 = 1 

where S 1 = {pi,Pi : Pi & {0, 1}, pi + p i — 1, i — 2, . . . , m — l} and the a^s are constants in 

[0,1]. 

Rewriting the right hand side of (3.2) gives 



7T7— 1 



7 = 1 



„ ,Pi_l«i-l+Piai \ I™ |P m _i«m-l+am 

•^"). I I 1^771 



7T7— 1 



Mm.O i m,ll^ C l ?7 - 2 / j I / Fll F2] M 
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x exp f N^ |a^| 2 it^ J du dx 

Integrating with respect to x gives 

n III _1_ 

\T n T n I <T n r n mS ^r 1 V^ / ,,-Hd-Ha 1 ^-Hd-Hp 2 a2TT^~Hd-H{p i _ 1 a i ^ 1 +p i a l ) 

\ 1 m,0 ~ 1 m,l\ ^ C 2 n Z^ / ! 2 11 « 

5 ^° m 7=3 

xn- / "- H( ^- iam - 1+am) ^. 

Assume that 

1 — i/d — Ha\ > 0, 1 — Hd — Hp 2 a2 > 0, 1 — Hd — H(p m _ 1 a m _i + a m ) > 

and 

1 — Hd — i7(pj_ 1 a i _i + PiOi) > for z = 3, . . . , m — 1. 

Then 

m 

Hd—i m(l—Hd)—H 5~1 ck* 

IC,o-Cil<c 3 n m ^(n6-na) & . 

Fix e > 0. We choose a± = 1 ~^ d — e if 1 — Hd < H. Otherwise, we let ot\ = 1. For 
i — 2, . . . , m, we choose a* = 1 ~^ d — e. With these choices of a^s, we obtain 

Hd-1 M TTJ , TT ^ i^+mHe if 1 - H d < H ; 

m h 77i(l — Hd) — H y a* — < , „, „„ 

2 l ; ^ 1 i-Hd-2g + ( m _ !)ff e m-Hd>H. 



Thus we can choose e such that m Hd 2 1 + m(l — iifcZ) — if J2 an = —7, and 



Hd-1 

T fioyi — siu) — si 

8=1 

, 1-ifd 



If _ r" I < r , „— T/'ft _ n \ m — 2 7 

which is the desired estimation. 

Step 2: Suppose now that k = 2. By the definition of I% ll and 7 mj2 > l-^mi — -^,2! is l ess than 
a constant multiple of 

I I \K-X2)\\Kx2-x-i)-T{x 2 )\ \\f{xi-x i+ i)\ exp ( — ^ V|xi| 2 u^)(iMda;. 



n 2 

i=3 i=l 

Using similar arguments as in Step 1, 

m—l 



m. 



Ci-C 2 l<c 5 n mi ^ i V / / \x 2 nx 3 r(~ [N^|x i+1 p<)| 

m 
/ ^// \ ^ I i9 2hf\ 

x exp I — y^ \Xi\ u i J du ax 

n m — / 
Sl Jo,.. 

m—l 
(' ' \ u .\~ Hd ^ H (P 1 -i a i~i+Pi a i)\\ u \-Hd-H(p m _ 1 a m -i+am) ^ u 



i=\ 



I I— Hd\ . \-Hd-Ha\\ \- H d- H (a2+P2,az) 
I u l I I ""21 I ""31 



m—l 

j=4 



where S\ = {pi, p i : Pi E {0, l},Pi + p i — 1, i — 3, ■ ■ ■ , m — l}. Then we can conclude as in 
Step 1. 

Step 3: Suppose that k is odd and 3 < k < m. Since k is odd, \I^ k _ ± — I™ k \ is less than a 
constant multiple of 

m 

n 1 ^ 



n 2 I [ \J{Xi-x i+1 

md I n V 



fc-1 
2 m 

I a l 

j=l - i=l 



x \f(xk - x k+1 ) - f(-x k +i)\( ] I l/O^j) I 2 ) exp ( - -^ S y^ j |a 

i=i 

Therefore, |/ mfc _! — / mfc | is less than a constant multiple of 

fc-i 

/> /» m 2 " ( 

^/ / (II lfe-^m)l)l^r fc (ni/K)| 2 )exp(-^^|x J |\^)^rfx. 

JW n J O m A — U 1 1 1' — 1 ; — 1 



fe-1 

2 m 



n m 



=fc+l i=l i=l 

Integrating with respect to x^s and u^s with i < k — 1 gives 

lmi-1 m.fcl 



^i — 3-i+lJI ) \ X k\ 



n n lib 

<c 7 {b-a)^-m n ( m -^)^ / ( ■ |/( Xi -x i+1 )|) 

J U (.m-k+l)i J 0mik V .^ 
m 

x exp I — y^ \Xi\ Mj ) dudx, 



where du = du k ■ ■ ■ du m , dx = dxk ■ ■ ■ dx m and 

O m ,k = {(u k ,...,u m ) :0<Ui< n(b- a),i = k,...,m}. 
Applying Step 1 and then doing some algebra, we can obtain 

ICk-i-Ckl^csn-^b-ar 1 ^^. 

Step 4 : The case when k is even and 4 < k < m is handled in a similar way. ■ 

4 Proof of Theorem 1.1 

The proof of Theorem 1.1 will be done in two steps. We first show tightness, and then estab- 
lish the convergence of moments. Tightness will be deduced from the following inequality. 

Proposition 4.1 For any < a < b <t and any integer m > 1, 

E [(F n (b) - F n (a)) 2m ] <C(b- a) m ^- Hd ^\ 
where C is a constant depending only on H , m, d and f . 

Proof. Note that E [(F n (b) — F n (a)) 2m ~\ = /^ m0 . Applying Proposition 3.1 repeatedly gives 

%nfl < ci n~\b - a)^- Hd ^ + Cl I^ 2m . (4.1) 

So it suffices to estimate l^mim- By the definition of l2m,2mi using the same notation as in 
the proof of Proposition 3.1, we obtain 



n n til 

(1 2m \ 

- -Var ( y^ Xj ■ (-B(sj) - B(si-i)) J J ds dx 

„ „ m 2m 

< c 3 n m ^ / / (J] \fM\ 2 ) exp ( - ^ J2 NV 2H ) dudx 
<c 4 {b-a) m{1 - Hd) ( / \f{x)\ 2 \x\~«dx) . (4.2) 

Combining (4.1) and (4.2) gives the desired result. ■ 

Next we shall prove the convergence of all finite dimensional distributions. That is, we 
shall prove that the moments of F n (t) converge to the corresponding ones of W(L t (0)). 

Fix a finite number of disjoint intervals (ai,bi\ with i = 1,...,N and bi < Oj+i. Let 

TV 

m = (mi, . . . , ttitv) be a fixed multi-index with m^ G N for i — 1, . . . , N. Set Yl m « = I 111 ! 



N 

and \~[ m i- = m - We need to consider the following sequence of random variables 

N 



G n = \l(F n (b l )-F n (a l )) r 



i=l 



and compute lim E (G n ). Note that the expectation of G n can be formulated as 

n— >oo 

N mi 



Iml(Hd-l) 

E(G n ) = m\n 2 






where 

fl m = {s G K |m| : na { < s\ < ■ ■ ■ < s\ TH < nb h 1 < i < N}. (4.3) 

Here and in the sequel we denote the coordinates of a point s E M' m ' as s = (s*), where 
1 < i < N and 1 < j < rrii. 

For simplicity of notation, we define 

J = {(i,j) :l <i < N,l < j <rm}. 

For any {i\,j\) and (22, J2) G ^cb we define the following dictionary ordering 

(hji) < («2, J2) 
if ix < i 2 or ^1 = 22 and j x < j 2 . For any (2, j) in J , under the above ordering, (i,j) is the 

C52 m k + j)-th element in J and we define #(z, j) = J^ m fc + j. 
fc=i fc=i 

Proposition 4.2 Suppose that at least one of the exponents m, is odd. Then 

lim E (G n ) = 0. 

Proof. Using Fourier transform, we see that E (G n ) is equal to 

<^k* n '^ LL (nn/fej)) e* P ( - ^(e» . B(s -))) **. 

Making the change of variables x^ — Yl yi ioi 1 < i < N and 1 < j < rrii, 

■ „ „ N mi 

E (G - } = 72^ nM ^ / ld / n n / w - 4 +1 ) 

(~ N mi v 

- 2 Var (EE s i ■ ( 5 ( s P - 5 (4-i))) ) rfs ^- 
i=i .7=1 ^ 



Applying Proposition 3.1, we obtain 



- -Var ( J2 E x l ' ( 5 ( s 5) - B ( S U))) ) rfscb ' 

8=1 j = l ' 



where J e = \{i,j) E Jo : Jf{i,j) is even} and 



/OmJ; if|m|isodd; 
1, if |m| is even. 

It is easy to see that lim E((7 n ) = when |m| is odd. We shall show lim E (G n ) = 

n— >oo n— >oo 

when I in I is even. In this case, 

N rrii 
X 



(1 N rrij v 

- 2 Var (EX>r ( 5 ( s P - 5 ( s 5-i))) <**<**■ 
i=l 1=1 ' 



Note that the right hand side of the above equality is positive. Using the local nondeter- 
minism property (2.2), 

I limE(G n )| ^dlimsupn^ 1 ^ 11 / I ( T\ \f(x))\ 2 ) 

/ N mi s. 

:= c\ lim sup I n . 

Assume that me is the first odd exponent. Integrating with respect to proper x^s and s*s 
gives 

I n <c 2 n H ^ sup / / / |/Vi +1 )| 2 (4-4-i)^ 

x exp ( - f\x{^(4 +1 - 4) 2 ") d < +l < dx[+\ 



Note that |/(x)| < c a \x\ a for a G [0, 1]. Choosing a e (^f^, ^f^) gives 
i n Sc 3 n sup / / / \x x I {s me - s me _ 1 ) 

x exp ( - «f \x^\\s^ - sl t r) dsf dsl ( dx^ 
< C4 n**" 1 sup / / (4+1 - sl t Y Hd ~ 2Ha (sl t ~ sl^Y™ A? 1 dsi le 

/nbi 
(nn </ \l-Hd-2Ha(J t \-Hd j J 

(na£ + i — s mt ) [s me - s mt _ 1 ) as mt 

L,-i 

where we used bi < a^+i in the last inequality. Therefore, 

I lim E(G n )| < c 6 lim n ^ Hd ~ 2Ha = o. 



Consider now the convergence of moments when all exponents m ; are even. 
Proposition 4.3 Suppose that all exponents m 8 are even. Then 



N 

|m| 



|m| / \ 

lim E(G„) = Cj d ||/||f I E( ] [ (W(JLfc(0)) - W(L a ,(0))) m M, (4.4) 

i=\ 

•where the expectation in the right-hand side of the above equation is given by formula (2.1). 
Proof. For any K > and £ = 1, . . . , |m|/2, we introduce the set 

D \m\,K = {( s l; • . . ,S| m |) e .D m : Si < S 2 < • • • < S| m |, S 2 i ~ S 2 £-l > K) , 

where D m is defined in (4.3). 

Taking into account the results proved in [3], the proof of the convergence (4.4) reduces 
to show that for all £ = 1, ... , |m|/2 

lim limsup w ^'^'^ Ef / \1 f(B( Si )) ds) = 0. (4.5) 

|m|,X 1 = 1 



In order to prove (4.5), set |m| = 2m. Then, it suffices to show that 

2m , 2m 

. limsuDn" 1 ^- 1 ) / / ( 



„ „ 2m 1 2m 

lim lim sup n m{Hd ' 1) / / ( TT /(a*)) exp f - -Var ( V a* ■ ^(s;))) ds dx = 0. 
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Using similar arguments as in the proof Proposition 3.1 we can write 

Mrnsup n m{Hd - 1) / ( T\ f( Xi )) exp ( - -Var (V x % ■ Bfa))) ds dx 

= limsnp^^)/ / (f[\f(z 2j )\ 2 ) 

n->00 JM?rnd J D i \ " / 

2m, K J— ^ 

(i 2m v 

- -Var ( ^ Zj • (-B(sj) - -B(sj-i))) J dscfe. 

The right hand side of the above equality is positive and less than or equal to 

„ „ m 2m 

limsupn^" 1 ) / / (]l\f(z v )\*)exp(-!%T l \z i \\8 i -8^ 1 ) VI )d8dz. 

2m, K J — -L 2=1 

Integrating with respect all ZiS and Si$ (i ^ 2£), the above limit is less than or equal to 



n /*oo 

ci / / \f(z 2e )\ 2 e-"f^ 2u2H dudz 
<c 2 / |^| 2 e-^^l 2 " 2H rfnrf^ 

1-Hd-2H 



= c 3 K 
This completes the proof since 1 — Hd — 2H < 0. ■ 

Proof of Theorem 1.1. This follows from Lemma 2.1, Propositions 4.1, 4.2 and 4.3 by 

the method of moments. 
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